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Spin properties of two interacting electrons in a quantum dot (QD) embedded in a nanowire with 
controlled aspect ratio and longitudinal magnetic fields are investigated by using a configuration 
interaction (CI) method and exact diagonalization (ED) techniques. The developed CI theory 
based on a three-dimensional (3D) parabolic model provides explicit formulations of the Coulomb 
matrix elements and allows for straightforward and efficient numerical implementation. Our studies 
reveal fruitful features of spin singlet-triplet transitions of two electrons confined in a nanowire 
quantum dot (NWQD), as a consequence of the competing effects of geometry-controlled kinetic 
energy quantization, the various Coulomb interactions, and spin Zeeman energies. The developed 
theory is further employed to study the spin phase diagram of two quantum- confined electrons 
in the regime of "cross over" dimensionality, from quasi-two-dimensional (disk-like) QDs to finite 
one-dimensional (rod-like) QDs. 

PACS numbers: 



I. INTRODUCTION 

Stimulated by recent success in coherent control of two- 
electron spin in laterally coupled quantum dots (QDs), 1 
the spin states of two interacting electrons in semicon- 
ductor QDs have received increasingly considerable at- 
tention. Accessible and engineerable spin states of few 
electrons in QDs thus have become one of the basic 
features required by the quantum information applica- 
tions in which electron spins are utilized as quantum 
bit. 2,3 For two-dimensional (2D) epitaxial QDs, magnetic 
field induced spin singlet-triplet (ST) transitions of two- 
electron ground states have been studied extensively for 
years4^&£ The underlying physics of the ST transitions 
is usually associated with the energetic competition be- 
tween quantized kinetic energies, the coulomb interac- 
tions, and spin Zeeman energies. Reversely switching 
the singlet and triplet spin states of a lateral two-electron 
QD is feasible by utilizing electrical control^ Moreover, 
it has been both theoretically and experimentally shown 
that more complex oscillating spin phases can be gen- 
erated either by reducing the lateral confinement or by 
increasing an applied magnetic field ^ 10 ! 11 ! 12 

Recently, the local-gate electrical depletio n 13 ! 14 ! 15 and 
the bottom-up grown technique s 16 ! 17 have been devel- 
oped for the fabrication of few-electron QDs embedded 
in a nanowire. These experimental developments open up 
an opportunity of exploring the cross over mechanisms 
from the 2D (disk-like) to the finite ID (rod-like) QD 
regimes. Such nanowire quantum dots (NWQDs) are ad- 
vantageous for geometrical control over a wide rage of as- 
pect ratio a (typically from a ~ 10 _1 to a > 1) j 16 t 17 The 
excellent versatility of shape and dimensionality makes 
NWQDs a suitable nanomaterial for scalable quantum 
electronics. Very recently, successful fabrication of sin- 
gle electron transistors made of InAs based gate-defined 
NWQDs and observations of the singlet-triplet transi- 



tions of two electrons in the QDs have been demon- 
strated. 18 How the highly tunable longitudinal confine- 
ment of NWQD affects and can be utilized to tailor the 
spin properties of few electrons in NWQDs are interesting 
subjects worth studying. 

The above experimental efforts motivate us to perform 
a theoretical investigation of the spin states of two elec- 
trons in InAs-based NWQDs— 8 by using a developed con- 
figuration interaction (CI) theory and exact diagonaliza- 
tion techniques^ The developed CI theory is based on 
the 3D parabolic model with arbitrary transverse and 
longitudinal confinement strength s 20 ! 21 and provides ex- 
plicit generalized formulations of the Coulomb matrix, 
and thus allows for straightforward and efficient numeri- 
cal or even semi-analytical implementation widely appli- 
cable for various cylindrically symmetric QDs. Our exact 
diagonalization studies of two-electron charged NWQDs 
with controlled geometric aspect ratios and longitudinal 
magnetic fields reveal fruitful features of spin singlet- 
triplet transitions, as a consequence of the competing ef- 
fects of geometry-engineered kinetic energy quantization, 
the various Coulomb interactions, and spin Zeeman ener- 
gies. The developed theory is further employed to study 
the spin phase diagram of two quantum-confined elec- 
trons in the regime of "cross over" dimensionality from 
quasi-2D (disk-like) QDs to finite ID (rod-like) QDs. 

This article is organized as follows: Section II describes 
the theoretical model and the developed configuration 
interaction theory for few-electron problems of three- 
dimensionally confining quantum dots. In Sec. Ill, we 
present and discuss the calculated results of magneto- 
energy spectrum, the ST transitions and geometry- 
engineered spin phase diagrams of two-electron charged 
quantum dots embedded in nanowires. Concluding re- 
marks are presented in Sec. IV. 
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where n, ra, g = 0, 1,2--- denote oscillator quantum 
numbers, Ez = g*\±BBs z is the spin Zeeman energy, 
oo± = uoh ±uo c /2 are in terms of the hybridized frequency 
uo h = (c^q +^/4) 1 / 2 and the cyclotron frequency uj c = 
eB/m*. The corresponding eigenstate |n, m, g) possesses 
the orbital angular momentum projection i z = h(n — m) 
and the parity P = 1 (P = —1) with respect to z-axis 
for an even (odd) q number. The wave function of the 
lowest orbital is given by 
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FIG. 1: (Color online) Single-electron energy spectrum as a 
function of aspect ratio a of a NWQD with fixed lateral con- 
finement huo = 13.3 meV at zero magnetic field obtained 
from the 3D parabolic model. The considered lateral con- 
finement strength huo — 13.3 meV corresponds to the cross 
section diameter Lo ~ 65 nm for a cylindrical InAs nanowire. 
The low- lying s-, p ± -, and p°-orbitals are relevant to a two- 
electron problem. The energy quantization for a short (long) 
NWQD with a < 1 ( a > 1) is characterized by the energy 
difference between the lowest and first excited orbitals huo 

{huj z ). 



II. MODEL 



with the characteristic lengths of the wave function ex- 
tents 1^ = ^h/2m*ujh and l z = y/h/2m*uj z . The 
wave functions of other excited states can be generated 
by successively applying the following defined raising 
operators^ 
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A. Single-particle model 



onto the ground state |0,0,0), i.e. 



We begin with the problem of a single electron in 
a NWQD with a uniform longitudinal magnetic field 
B = (0,0,5), which is described by the single-electron 
Hamiltonian, 
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-(p + eAy + V c (x,y,z)+H z . 
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Here the first term indicates the term of kinetic energy 
with A = (B/2)(—y,x,0) being the vector potential in 
symmetric gauge, m* the effective mass of electron and e 
the charge of an elctron. The second term is the confining 
potential of NWQD modeled by 



1 



V c (x, y, z) = ^m* [cjq (x 2 + y 2 ) + u 2 z z 2 ^ 



(2) 



with uq and uo z parametrizing, respectively, the trans- 
verse and the longitudinal confining strength. The last 
term is the spin Zeeman energy Hz = g* HbBs Zi in terms 
of the z-component of electron spin s z = ±1/2, the effec- 
tive Lande ^-factor of electron g* and the Bohr magneton 
He- The single-particle Hamiltonian (pQ) leads to the ex- 
tended Fock-Darwin single-particle spectrum 
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The diameter of cross section of bottom-up synthe- 
sized nanowire is typically ~ 50 — 70 nm. By contrast, 
the length of a QD in a nanowire, defined by imposed 
electrodes or heterostructure potential barriers, is highly 
tunable over a wide range from 10 to 300 nm^ For char- 
acterizing the geometry of a NWQD, it is convenient to 
define the parameter of aspect ratio, 
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according to the characteristic length of the lowest or- 
bital wave function based on the 3D parabolic model. A 
rod-like (disk-like) NWQD is characterized by the value 
of aspect ratio a > 1 (a < 1), where the longitudinal ex- 
tent of the electron wave function is longer (shorter) than 
the transverse one on the cross section of the nanowire. 
Notably, the effective aspect ratio a = I z /Iq defined here 
is not but very close to the geometric aspect ratio a geom , 
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(3) section diameter (length) of NWQD. 



L^/Lo with Lq (L^) being the cross 
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FIG. 2: (Color online) Schematic illustration of the electronic structures, consisting of few relevant low lying orbitals (one s- 
and three p-orbitals), of long rod- like NWQDs [(a)(c)(e)] and short disk- like NWQDs [(b)(d)(f)] with or without longitudinal 
magnetic field B and including or excluding the spin Zeeman splitting Ez (g* = or g* ^ 0). (a) a > 1 and 5 = 0; (b) a < 1 
and B = 0; (c) a > 1, B ^ and = 0; (d)a < 1, B ^ and g* = 0; (e)a > 1, 5 ^ and g* ^ 0; (f)a < 1, B ^ and g* ^ 0. 



Figure [T] presents the calculated single-electron energy 
spectrum as a function of aspect ratio a for a N WQD 
with fixed lateral confinement hcoo = 13.3 meV at zero 
magnetic field according to Eq.([3]). The chosen param- 
eter of lateral confinement hcoo = 13.3 meV is deter- 
mined by fitting the numerically calculated energy sepa- 
ration between the two lowest single-electron orbitals of 
a cylindrical InAs/InP NWQD of cross section diame- 
ter Lo = 65 nm by 3D finite difference simulation. In the 
simulation, the Schrodinger equation for a single electron 
confined in a 3D cylindrical potential well is solved by us- 
ing finite difference method, with the used parameters: 
the effective mass m* = 0.023rao of electron for InAs and 
the InAs/InP band edge offset Vb = 0.6 eV as the barrier 
height of the confining potential ] 17 i 22 

In a two-electron (2e) problem, the most relevant or- 
bitals are the two lowest ones because the kinetic energy 
difference between the two orbitals is the main energy 
cost, in competition with the coulomb or spin Zeeman 
energies, for a spin triplet state to be the ground state 
of two-electron. By convention, we from now on name 
the lowest single electron state |n, m, q) = |0,0,0) as s- 
orbital, and the next three p-shell states |0,0,1), |1,0,0), 
and |0, 1, 0) asp -, and p~ -orbitals, respectively. Ac- 
cording to Eq. (j3]), the energy of the lowest s-orbital is 
explicitly given by 

e s , Sz = ~ + ficj- + Tujz) + g* ^bBs Zj (8) 



and those of the three p-shell orbitals are respectively 
given by 

e p°,s z = e s,s z + fi^z, 
e p+,s z = e s,s z + 

e P -,s z = £s,s z + ftu— (9) 

For B = 0, we have e SySz = huo (l + l/2a 2 ), e p o^ Sz = 
e s ^ Sz + hcuo/a 2 , and e p + ySz = e p -^ Sz = e s , Sz + fkj accord- 
ing to Eqs.(0 and ©. Here, the p + - and p~ -orbitals 
are degenerate with the same energy separation from the 
s-orbital, hw± = Hujq, while the p°-orbital is energeti- 
cally higher than s-orbital by Hlj z = hcoo/a 2 . Obviously, 
p° (p ± ) is the second lowest orbital for a long (short) 
NWQD with a > 1 (a < 1) at zero magnetic field, as 
shown in Fig. [TJ For a symmetric NWQD with a = 1, 
the p°- and p ± -orbitals form a 3-fold orbital-degenerate 
shell. Figure [2] (a) [(b)] schematically depicts the low- 
lying orbitals of a long [short] NWQD with a > 1 [a < 1] 
at zero magnetic field. 

Applying a longitudinal magnetic field onto a cylindri- 
cal NWQD breaks the degeneracy of p + - and p~-orbitals. 
The orbital Zeeman effect lowers (raises) the energy level 
of the p~(p + )-orbital from huo to huj- (Hlj+). Thus, if 
a long NWQD is subjected to a sufficiently strong mag- 
netic field, the second lowest orbital of the dot could be 
changed from the p° to p~ . By contrast, the second low- 
est orbital of a short NWQD is always the p~ -orbital. 
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Therefore, the characteristic energy quantization of the 

1/2 

p~-orbitals, Hlj- = h [cOq + o; 2 /4) — hjj c /2, is often a 
key parameter for a short NWQD or a moderately long 
NWQD with strong magnetic field. Considering wide- 
band gap materials such as GaAs, the ^-factors are usu- 
ally small and the spin Zeeman effect on the energy shift 
of orbital is negligible. Figure [2](c) [(d)] depicts the in- 
dependent electronic orbitals of a long [short] NWQDs, 
where vanishing spin Zeeman splitting is assumed (g* = 
is set). 

For a low energy gap material with larger g* , like InAs, 
the spin Zeeman effect could be significant in the spin ST 
transition of two-electron QD. Figure EJe) [(f)] schemati- 
cally shows the spin-resolved electronic orbitals of a long 
[short] NWQDs with B ^ and g* ^ by the spin 
Zeeman splitting 2Ez- With the spin Zeeman effect, all 
the spin-up (spin-down) orbitals are energetically lowered 
(raised) by E z = g*HBB/2 according to Eq. ([3]). If the 
applied magnetic field or the ^-factor of material is so 
large that the spin Zeeman splittings exceed the kinetic 
energy quantization of QD, both of the two lowest single- 
electron states are the spin-up ones and the ground state 
of the 2e dot is ensured to be a spin triplet state simply 
according to spin Pauli exclusion principle. 

In this work, the following formulation for the ^-factor 
of an InAs-based QD is adopte d 17 ! 23 



3 Ef {Ef + A so ) 



(10) 



where Ef is the effective energy gap of semiconductor 
QD, g = 2.0 is the Lande ^-factor for free electron, Aso 
is the spin-orbit splitting in the valence band, and P is 
the parameter of interband transition matrix element. 23 
Here, the effective energy gap of a QD can be esti- 
mated as Ef = £kuik + 6s ^ ^ where ^buik is the bulk 

energy gap and e SjSz is the quantization energy of the 



lowest electronic orbital of the QD with B = measured 
from the conduction band edge. For InAs-based QDs, we 
take the following parameter values: E^ ulk = 460 meV, 
A S o = 390 meV, P 2 = 21.5 eV. 17 Accordingly, the value 
of g* for a symmetric NWQD with Lq = L 2 = 65 nm is 
estimated as large as g* « —11^ 



B. Interacting few-electron model 

The interacting Hamiltonian of few electrons in a 
NWQD can be expressed in the form of second quan- 
tization as 

H = Y] £iac\(jC-io- 

+ \ E (ij\V\kl)clcl,c k(T ,c la , (11) 

ijkl,aa' 

where z, j, /c, I denote the composite indices of single elec- 
tron orbitals such as \i) = \rii,mi,qi), c\ a (q ct ) the elec- 
tron creation (annihilation) operators, and a = ± the 
electron spins s z = ±|. The first (second) term on the 
right hand side of Eq. ([TT]) represents the kinetic energy 
of electrons (the Coulomb interactions between electrons) 
and the Coulomb matrix elements are defined as 

(ij\V\kl) = -^-J J dndr 2 ^(ri)^(r 2 ) 

x, r^fe(r 2 M(ri), (12) 

|ri - r 2 | 

where k is the dielectric constant of dot material. For 
InAs material, we take k = 15.15. After lengthy deriva- 
tion, one can obtain the generalized Coulomb matrix el- 
ements for the case of a > 1: 



{nimiqi\njmjqj\V\nkm k qk\nimiqi) = 
1 N \ S R , 
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2 r ( 3+2W 2 iFijl + u, ; ;l-x) , (13) 



where we define u = m i +mj+ni+n k -(p 1 +p 2 +P4,+P5), (ni+nj) = -(4,i + Rr = ( m i +m k ) - (n ( +n k ) = 
v = (qi+qi + qj + q k ) - 2(p 3 + p 6 ), R L = (rrn + mj) - 
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— (£z,i +^,fe)? x = ^z/ujh, and 2^1 is the hypergeometric 
function. The ^-functions 5 qi + q .+ qi + qkieven and S RljRr 
in the formulation ensure the conservation of the parity 
with respect to z-axis and the ^-component of angular 
momentum of system L z , respectively. The formulation 
of Eq. ([T3j) is confirmed by computing the Coulomb in- 
tegral numerically. 

For short NWQDs with a < 1, the formulations of the 
Coulomb matrix elements are obtained by simply taking 
Euler's hypergeometric transformation for the hypergeo- 
metric function in Eq. ([T3]h i.e., replacing 
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The generalized formulations for the Coulomb matrix 
elements based on the 3D asymmetric parabolic model 
are probably for the first time derived, which allows for 
straightforward implementation of the CI theory and is 
widely applicable to arbitrary 3D confining semiconduc- 
tor nanostructures. 



C. Exact diagonalization 

Based on the CI theory presented above, we follow the 
standard numerical exact diagonalization procedure to 
calculate the energy spectrum of N e interacting electrons 
in a NWQD . 19 The numerically exact results are ob- 
tained by increasing the numbers of chosen single electron 
orbital basis and the corresponding 7V e -electron configu- 
rations until a numerical convergence is achieved. In the 
full configuration interaction (FCI) calculation for a 2e 
problem, we usually take the number of single electron 
orbitals typically from 20 to 26 and that of the corre- 
sponding 2e configurations from 190 to 325 to have a 
satisfactory numerical convergence. 



III. NUMERICAL RESULTS AND DISCUSSION 

A. Magnetic-field induced ST transitions 

Let us first consider two interacting electrons in a rod- 
like NWQD with the aspect ratio a = 3 and the trans- 
verse confining strength hcoo = 13.3 meV using FCI cal- 
culation. The low-lying magneto-energy spectrum of the 
two-electron NWQD is shown in Fig. [3l which consists of 
a spin singlet state branch, labeled by S, and three triplet 
state branches split by the spin Zeeman energy, labeled 
by T^q, T lo and T^ according to the z-component of 
total spin (S z = +1, S z = and S z = —1), respec- 
tively^^ Since usually only triplet states with S z = +1 
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are involved in ST transitions, we shall use T L |l z | to 
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FIG. 3: (Color online) Magneto-energy spectrum of two in- 
teracting electrons in a NWQD with transverse confining 
strength hu>o = 13.3 meV and aspect ratio a = 3. 



denote the triplet states with angular momentum L z 
through out this article, skipping the superscript + of 
T^| L | for brevity. 

The main configurations of the two-electron ground 
states around the critical magnetic field are schematically 
shown in the lower right corner of Fig. [3j In the weak 
magnetic field regime (B < £>st lo ~ 0.9 T), the two elec- 
trons in the NWQD mainly occupy the lowest s-orbital 
simply following the Aufbau principle, and form a spin 
singlet ground state. With increasing £?, the triplet state 
Tlo is more energetically favorable than the singlet state 
because of the increasing spin Zeeman energy, the re- 
duced Coulomb repulsion, and exchange energy between 
the two spin polarized electrons. A crossing of the singlet 
branch and the triplet state branch Tlo is observed at the 
critical magnetic field £?st lo = 0.9 T. Such magnetic-field 
induced ST transitions are attributed to the energetic 
competition between single particle energy quantization, 
the spin Zeeman energy, and the various Coulomb inter- 
actions including the direct, exchange, and correlation 
interactions as well. 10 

Other weak spin-related terms, such as the spin-orbital 
coupling (SOC) with 1-2 order of magnitude smaller 
than the kinetic quantization of QD are neglected in the 
Hamiltonian of Eq. ([TT]) . The SOC mixes the spin of the S 
and Tlo states and creates an anti-crossing of the S- and 
TLo-branches around the Bst with a small energy gap, 
typically only ~ 0.1 — 0.5 meV as observed in previous 
experiments^ 



B. Spin phase diagram 

Figure [H shows the calculated spin phase diagrams of 
the two-electron ground state of the NWDQs with a fixed 
cross section diameter (fixed hcoo = 13.3 meV) but var- 
ious lengths (various hu z ) with respect to the applied 
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a 

FIG. 4: (Color online) Spin phase diagrams of two-electron 
NWQDs of lateral confinement hu>o — 13.3 meV with re- 
spect to tunable magnetic field B and aspect ratio a. The 
phases are distinguished by the curves of critical magnetic 
field Sst obtained from non- interacting (black dotted), PCI 
(blue dashed), and FCI (red solid) calculations. 



magnetic field B and the aspect ratio a. Three phases 
(S, Tlo, and Tli) are distinguished by the curves of crit- 
ical magnetic field £>st in Fig. [H Correspondingly, the 
main configurations of the 2e ground states are depicted 
inside the colored regions of the phases. To identify the 
various underlying mechanisms in the phase diagrams, 
including the spin Zeeman effect and the inter-particle 
Coulomb interactions, the spin phase diagrams are cal- 
culated by using non-interacting, full CI, and partial CI 
calculations, respectively. 

In the non-interacting calculation, the coulomb inter- 
actions are artificially disabled and the considered ST 
transitions are induced only by the spin Zeeman effect. 
The comparison between the results of non-interacting 
and FCI calculations allows us to distinguish effects of the 
Coulomb interaction and spin Zeeman coupling on the ST 
transitions. In particular, to highlight the Coulomb cor- 
relation effect, a partial configuration interaction (PCI) 
calculation is also performed for the spin phase diagrams, 
in which only the lowest energy configuration is taken as 
the sole basis and the couplings from higher energy con- 
figurations are excluded. 

The essential features of the phase diagrams can be 
realized based on the non-interacting picture. For a not 
very long (small or moderate a) NWQD with weak 5, 
the 2e ground state is likely to be the spin singlet state 
S, simply following Aufbau principle (the yellow region 
in Fig. Starting from the singlet phase S, the two- 
electron ground state of a NWQD might be switched to 
the spin triplet phases (the pink region Tlo or the cyan 
region Tli) by increasing either a or B (see the horizontal 
and vertical dashed lines with arrows in Fig. [H for the 
guidance of eyes). 

Following the horizontal dashed line, the longitudinal 



energy quantization hco z is decreased by the increase of a. 
With the addition of spin Zeeman term, the spin-up level 
of p°-orbital could become even lower than the spin-down 
level of s-orbital if the decreasing hw z is so small as that 
hw z < 2\E Z \ (see the difference between the schematic 
configurations for the S and Tlo states). In this situa- 
tion, the 2e ground state can transit to the spin triple 
states Tlo, simply following spin Pauli exclusion princi- 
ple. On the other hand, the transition of a 2e ground 
state of NWQD from the singlet state S to the triple 
one Tli is shown also possible by increasing the strength 
of applied magnetic field. Following the vertical dashed 
line, increasing B reduce the energy separation between 
the s- and ^"-orbital levels, i.e. hco-. Similar to the 
case of S-Tlo transition, a S-Tli transition can happen 
as the decreased Hlu- is so small as that hw- <2\E Z \ . In 
the non-interacting picture, £?st lo * s ex Pli c itly given by 
^st lo = ^o/^*Mb^ 2 , showing a quadratic decay with a, 
while the critical magnetic field £?st L i ^ or S-Tli transi- 
tions is dependent only on Hluo and remains nearly con- 
stant in the B — a plot. 

The Coulomb interactions are shown to reduce the 
singlet phase area in the diagrams from the compari- 
son between the non- interacting and CI results. For 
example, the segment of vertical solid line at a = 3 in 
Fig. 2] indicates that the critical magnetic field is signif- 
icantly reduced from £?sT L0 (Non — interacting) = 3.2 T 
to £>st L0 (FCI) = 0.9 T as the Coulomb interactions are 
taken into account. This is because the spin triplet states 
gain additional negative exchange energies while the sin- 
glet state does not. We also notice that the £?st L i f° r ^ ne 
S-Tli transition no longer remains constant but slightly 
increases with increasing a because the strength of the 
coulomb interactions is reduced by the increase of dot 
volume. 

Basically, the results obtained from the FCI and PCI 
calculations have similar features except for those in the 
regime of high a (a > 3). While the PCI calculation 
shows the vanishing £?st lo f° r a ^ 3, the FCI calculation 
yields the always non-zero £?st lo - This means that the 
Coulomb correlations energetically favor the spin singlet 
state as ground state and become more pronounced in 
long NWQDs. 



C. Crossover from disk-like to rod-like QDs 

The spin phase diagrams of Fig. [4] suggest that pur- 
posely accessing a specific spin phase of two-electron is 
feasible through the geometrical control of NWQDs. For 
instance, the ground state of a two-electron NWQD can 
be switched from the singlet S to the triplet state Tlo by 
increasing the aspect ratio a at the fixed B = 5 T (trace 
the horizontal dashed line in Fig. [4]). 

Figure [5] presents the spin phase diagrams of two- 
electron NWQDs with respect to the lateral and lon- 
gitudinal confinements, parametrized by Hluq and hw Zl 
respectively, in a fixed magnetic field B = 5 T for (a) 
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FIG. 5: (Color online) Spin phase diagrams of doubly charged 
NWQDs with respect to the lateral and longitudinal confine- 
ments, parametrized by Huoq and hcu z , respectively, in a fixed 
magnetic field B — 5 T for (a) non-interacting two electrons 
with g* 7^ 0, (b) interacting two electrons with g* ^ and 
(c) interacting two electrons with g* = 0. 



( a ) a > 1, g # 

s a zK 

P — — h® + 

S -A- 


fl > 1, g* = ( c ) 

1 S 

P i A n& + 


Tlo ~ 
P 

S 5 


Tlo 

P 


Tli 

~*~ -A- 


Tli 



( b ) fl<l,g*5t0 

s z - 
p — 

S 

-A- 


a<l,g=0 ( d ) 
P 


Tli _ = 

d — 

" P is . 

-*- 


Tli 


f =d -n 

* p is - 




^ s ^ 


f - d "IP -s - 


Cl ^- p s S L4 



FIG. 6: (Color online) Relevant two-electron configurations 
possibly being the main components in the ground states of 
NWQDs in an uniform magnetic field B — 5 T for (a) a > 1 
and g* ^ 0, (b) a < 1 and g* ^ 0, (c) a > 1 and g* = 0, and 
(d) a < 1 and = 0. 



non-interacting two electrons with g* ^ 0, (b) interact- 
ing two electrons with g* / 0, and (c) interacting two 
electrons with g* = 0. In Figs. EH(a) and (b), we present 
the relevant two-electron configurations to the spin phase 
diagrams of Figs. 03(a) and (b) with the inclusion of spin 
Zeeman effect (g* ^ 0), while in Figs. EJc) and (d) we 
present the relevant two-electron configurations to the 
spin phase diagrams of Fig. [5](c) for g* = 

The non-interacting spin phase diagram is first shown 
in Fig. (5](a) in order to identify the spin Zeeman effect 
and also contrast the Coulomb interaction effects on the 
interacting spin phase diagrams presented in Fig. [51(b). In 
the non-interacting case, the features of the spin phases 
of Fig. [5](a) are purely determined by the competition 
between geometry-dependent quantized electronic struc- 
tures of dots and the spin Zeeman splitting, which is 
nearly a constant here created by the fixed B. Three 
distinctive spin phases, S, Tlo, and Tli, are marked in 
different colors in Fig. [51(a). In the yellow region where 
both Hujo {huj-) and hw z are large, the kinetic quanti- 
zations in both longitudinal and lateral directions are 
stronger than the spin Zeeman splitting and S remains as 
a ground state. Reducing the longitudinal confinement, 
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fouj z , can lead to the S-Tlo (from the yellow to the pink 
region) transition as ftw z < 2\E^\. Similarly, reducing 
the transverse confinement leads to the S-Tli transition 
as hw- < 2^z(from the yellow to the light cyan region). 

Compared with Fig. Ufa), the interacting spin phase 
diagram of Fig. [51(b) shows the following additional fea- 
tures: 

(i) Larger areas of both Tlo and Tli phases are ob- 
served because of the additional negative exchange ener- 
gies and the reduced direct Coulomb repulsions gained 
by the triplet states. 

(ii) A NWQD with hujQ ^12 meV could experience 
a three-phase transitions from Tli (cyan) to S (yellow), 
and then to Tlo (pink) with increasing the length of wire, 
from hw z > 25 meV to hw z < 5 meV (see the vertical line 
positioned at hcoo = 12 meV in Fig. 0(b)). 

(iii) In the regime of small huoo and large hw z (i.e. flat 
quasi-2D dots with a <C 1), a series of transitions from 
the spin single states to various triplet states, Tli, Tl3, 
Tls, etc. [see Fig. EJb)] and a staircase increase of total 
orbital angular momentum are observed with reducing 
the lateral confinement hcoo. 

In the weak laterally confining regime, few electrons in 
the quasi-2D QD in a high magnetic field successively fill 
the orbitals with negative ^-projection of orbital angu- 
lar momentum, i.e. the orbitals of lowest Landau level 
(LLL), with small kinetic energy separation hcj-. The 
inter-particle Coulomb interactions thus become partic- 
ularly pronounced among the particles on the nearly de- 
generate LLL orbitals with alomost quenched kinetic en- 
ergies. In order to minimize the coulomb repulsion, the 
particles on the quasi-degenerate orbitals tend to spread 
the occupancy of orbitals as far as possible, but in com- 
petition with the cost of increase of kinetic energy. As a 
result, with reducing Hujo or increasing £?, the total an- 
gular momentum of two-electron increases, as previously 
discussed by Wagner et al£ for gated 2D QDs. 

Figure [5tc) shows the phase diagram of two interacting 
electrons calculated by FCI method but with the vanish- 
ing spin Zeeman term, i.e. g* = 0. This allows us to 
distinguish the effects of spin Zeeman energy and the 
coulomb interactions on the spin phase diagram of Fig. 
03(b), and also to study the spin phases of QD made of 
a material with small g* such as GaAs. Without spin 
Zeeman splitting, the significant features of Fig. [5](c) 
are completely determined by the many-body effects and 
geometry-engineered electronic structures of NWQDs. 

In the a > 1 regime, unlike the result shown in Fig. 
E](b), the Tlo phase disappears and naturally there is no 
S-Tlo transition observed. This is because the Coulomb 
correlations that energetically favor spin-singlet state as 
mentioned previously, become dominant and compen- 
sate the negative exchange energy gained by the Tlo 



states However, in the small huoo regime, an addi- 
tional singlet-triplet state oscillation with decreasing huoo 
is observed. Compared with Figure 03(b), the difference 
is the emergences of various singlet states between the 
triplet phases. This is due to the removal of spin Zee- 
man splittings, which energetically favor only the triplet 
states. Such a singlet-triplet state oscillation is evidenced 
as a main feature of a flat 2D QD with small spin Zeeman 
effect, as shown both theoretically 9 and experimentally 12 
in the previous studies. 



IV. SUMMARY 

In conclusion, we present exact diagonalization stud- 
ies of spin phase transitions of two electrons confined in 
nanowire quantum dots with highly tunable aspect ra- 
tio and external magnetic field. A configuration inter- 
action theory based on a 3D parabolic model for such 
three dimensionally confining QDs is developed, which 
provides generalized explicit formulation of the Coulomb 
matrix elements and allows for straightforward imple- 
mentation of direct diagonalization. The exact diago- 
nalization study reveals fruitful features of spin ST tran- 
sitions with respect to the tunable geometric aspect ratio 
and applied magnetic field. 

For disk-like QDs, the ST transition behaviors may be 
dominated by the spin Zeeman, the direct- Coulomb, and 
the exchange energies. The pronounced Coulomb corre- 
lations are identified in rod-like QDs with aspect ratio 
a > 3, which energetically favor singlet spin states and 
yield the always non-zero critical magnetic fields of ST 
transitions. The developed theory is further employed to 
study spin phase diagram in the dimensional "cross over" 
regime from the 2D (disk-like) QDs to finite ID (rod-like) 
QDs. In the 2D disk-like QD regime, various distinctive 
spin phases are emerged under the conditions of appro- 
priate lateral confinement strength and magnetic fields. 
In the rod-like QD regime, switching the ST transitions 
is shown feasible by controlling both lateral and/or lon- 
gitudinal confinement strength. 



V. ACKNOWLEDGEMENT 

This work was financially supported by the National 
Science Council in Taiwan through Contracts No. NSC- 
98-2112-M-009-011-MY2 (SJC) and No. NSC97-2112- 
M-239-003-MY3 (CST). The authors are grateful to the 
facilities supported by the National Center of Theoretical 
Sciences in Hsinchu and the National Center for High- 
Performance Computing in Taiwan. 



* Electronic address: sjcheng@mail.nctu.edu.tw 



1 J. R. Petta, A.C. Johnson, J. M. Taylor, E. A. Laird, A. 



9 



Yacoby, M. D. Lukin, CM. Marcus, M. P. Hanson, and 
A. C. Gossard, Science 309, 2180 (2005). 

2 D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

3 J. Fischer, M. Trif, W. A. Coish, and D. Loss, Solid State 
Comm. 149, 1443 (2009). 

4 L. P. Kouwenhoven, T. H. Oosterkamp, M. W. S. Danoe- 
sastro, M. Eto, D. G. Austing, T. Honda, and S. Tarucha, 
Science 278, 1788 (1997). 

5 L. P. Kouwenhoven, D. G. Austing, and S. Tarucha, Rep. 
Prog. Phys. 64, 701 (2001). 

6 S. M. Reimann, and M. Manninen, Rev. Mod. Phys. 74, 
1283 (2002). 

7 C. Ellenberger, T. Ihn, C. Yannouleas, U. Landman, K. 
Ensslin, D. Driscoll, and A. C. Gossard, Phys. Rev. Lett. 
96, 126806 (2006). 

8 J. Kyriakidis, M. Pioro-Ladriere, M. Ciorga, A. S. Sachra- 
jda and P. Hawrylak, Phys. Rev. B 66, 035320 (2002). 

9 M. Wagner, U. Merkt, and A. V. Chaplik, Phys. Rev. B 
45, 1951 (1992). 

10 P. Hawrylak, Phys. Rev. Lett. 71, 3347 (1993). 

11 B. Partoens, A. Matulis, and F. M. Peeters, Phys. Rev. B 
59, 1617 (1999). 

12 Y. Nishi, Y. Tokura, J. Gupta, G. Austing, and S. Tarucha, 
Phys. Rev. B 75, 121301 (R) (2007). 

13 C. Fasth, A. Fuhrer, M. T. Bjork, and L. Samuelson, Nano 
Lett. 5, 1487 (2005). 



A. Pfund, I. Shorubalko, R. Leturcq, and K. Ensslin, Appl. 
Phys. Lett. 89, 252106 (2006). 

A. Pfund, I. Shorubalko, K. Ensslin, and R. Leturcq, Phys. 
Rev. Lett. 99, 036801 (2007). 

M. T. Bjork, C. Thelander, A. E. Hansen, L. E. Jensen, 
M. W. Larsson, L. R. Wallenberg, and L. Samuelson, Nano 
Lett. 4, 1621 (2004). 

M. T. Bjork, A. Fuhrer, A. E. Hansen, M. W. Larsson, L. 
E. Froberg, and L. Samuelson, Phys. Rev. B 72, 201307(R) 
(2005). 

C. Fasth, A. Fuhrer, L. Samuelson, V. N. Golovach, and 

D. Loss, Phys. Rev. Lett. 98, 266801 (2007). 

A. Wensauer, M. Korkusihski, and P. Hawrylak, Solid 
State Comm. 130, 115 (2004). 

M. Dineykhan, and R. G. Nazmitdinov, Phys. Rev. B 55, 
13707 (1997). 

J. T. Lin, and T. F. Jiang, Phys. Rev. B 64, 195323 (2001). 
M. T. Bjork, B. J. Ohlsson, C. Thelander, A. I. Persson, 
K. Deppert, L. R. Wallenberg, and L. Samuelson, Appl. 
Phys. Lett. 81, 4458 (2002). 

C. Hermann and C. Weisbuch, Phys. Rev. B 15, 823 
(1977). 

R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha 
and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 
(2007). 



